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ABSTRACT 

Type I X-ray bursts and superbursts on neutron stars release sudden and intense radiation fields into their 
surroundings. Here, we consider the possible effects of these powerful explosions on the structure of the 
accretion disk. The goal is to account for the apparent evolution of the innermost regions of the accretion 
disk around 4U 1820-30 during a superburst. Three different processes are considered in detail: radiatively 
or thermally driven outflows, inflow due to Poynting-Robertson drag, and a structural change to the disk by 
X-ray heating. Radiatively driven winds with large column densities can be launched from the inner disk, 
but only for L/L-Edd ^ 1, which is expected to be obtained only at the onset of the burst. Furthermore, the 
predicted mass outflow rate is less than the accretion rate in 4U 1820-30. Estimates of the Poynting-Robertson 
or radiative drag timescale shows that it is a very efficient means of removing angular momentum from the gas. 
However, the analytical results are likely only applicable at the innermost edge of the disk. X-ray heating gives 
a change in the disk scale height that is correlated with the blackbody temperature, as seen in the evolution 
during the 4U 1820-30 superburst. If this change in the scale height can alter the surface density, then the 
viscous time (with a ~ 0.03-0.2) is the closest match to the 4U 1820-30 results. We expect, however, that 
all three processes are likely ongoing when an accretion disk is subject to a sudden heating event. Ultimately, 
a numerical simulation of a disk around a bursting neutron star is required to determine the exact response of 
the disk. Magnetic truncation of the accretion flow is also considered and applied to the 4U 1820-30 X-ray 
reflection results. 

Subject headings: accretion, accretion disks — radiation mechanisms: thermal — stars: individual (4U 1820- 
30) — stars: neutron — X-rays: bursts 



1. INTRODUCTION 

Accretion disks around neutron stars are often subject to sudden intense radiat i on fields due to Type I X-ray bursts 
jLewin. van Paradi is & Taam|ll993 Ul995HBildstenlll9 98: Strohmaver & Bildsten 2 0031 ICum ming 20 041), or , on rare occasions, 
superbursts (Cornelisse et al. 2000; Cumming & Bildsten 2001; Strohmaver & Brown 2002; Kuulkers 2004). These events can 
reach the Eddington luminosity, and thus deposit a large amount of energy in the accretion disk over a short amount of time 
(tens of seconds for Type I bursts; hours for superbursts). The accretion disk structure may be altered by this sudden injection of 
energy, as suggested by the evoluti on of X-ray reflection features observed from the low-mass X-ray binary (LMXB) 4U 1820- 
30 during a superburst (Ballantvne & Strohmaver 2004). In this source, modeling of Fe Ka line and edge features produced 
by the disk indicated that the reflecting part of the disk moved outward from ~ 20 r g (r g = GM/c 2 , where M is the mass of 
the neutron star) to ~ 100 r g in ~ 1000 s, before eventually returning to ~ 10 r g . This behavior was puzzling because the 
chan ge in the reflecting radius w as not correlated with the observed X-ray flux, but with the hardness of the blackbody spectrum. 
Ballantvne & Strohmaver ( 2004) therefore proposed a model in which the surface density of the inner accretion disk was altered 
by heating from the incident X-rays; however, alternative models have not yet been carefully considered. 

In this paper, we examine in detail a number of different physical processes which may alter the structure of accretion disks 
during a sudden heating event. The physics are treated as generally as possible, but quantitative comparisons are made with the 
4U 1820-30 results in order to test the individual models as possible explanations for the inferred behavior. 

We begin in the next section by reviewing the observed properties of 4U 1820-30. The paper is then broken into three sections 
based on the different possibilities for the accretion disk response. In §[3]we consider outflow mechanisms, such as thermally and 
radiatively driven winds. §|4]contains a study of inflow onto the central object via radiative drag. Changes to the disk structure 
that do not involve inflow or outflow are the subject of §|5] The results are gathered and discussed in §[6] Finally, conclusions are 
drawn in §0 

2. OVERVIEW OF 4U 1820-30 

4U 1820-30 is one of the most compact LMXBs know n with an orbital period of only 11.4 m inutes 
( Stella, Priedho rskv & Whitell987l). Ev olutionary models of the system ( Verb unlll9"87llBailvn & Grindlavll 9871: 1 Rappapo rt et alJ 
1987; Rasio. Pfahl & Rappaport 2000) all predict tha t the accreting material is primarily helium. Type I X-ray b ursts from 
4U 1820-30 have been observed for nea rly 30 years JGrindlav et alJll976l lHaberl et al.lll987l iKuulkers et al.l [2003). and have 
recurrence times between 2 and 4 hours. Cumming (2003) presents further details and models of the Type I X-ray bursts from 
4U 1820-30. 

Although the system resides within the metal rich globular cluster NGC 6624, the distance to 4U 1820-30 is still not precisely 
known. The difficulty arises because of uncertainties in the interstellar extinction along the line-of-sight, and in the metallicity 
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relation used to correct the absolute magnitude of horizontal branch stars. Thus, using the s ame technique. | K uulkers et alJ (120031) 
derived the distance to 4U 1820-30 of 7.6±0.4 kpc, in agreement with an earlier estimate (Rich. Minniti & Liebert 1993), while 
IVacca. Lewin & van Paradiis (1986) found it to be 6.4±0.6 kpc. Most recently, Shaposhnikov & Titarchuk (2004) modeled 
the light-curve and spectrum of a photospheric radius expansion (PRE) burst from 4U 1820-30 and determined a distance of 
~ 5.8 kpc. 

The Eddington luminosity at the photosphere of a neutron star with a helium atmosphere as determined by a distant observer is 



L Ed d = 2.5 x 10 
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Mm 
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(1) 



where i?» is the radius of the photosphere (iLewin et alJl993|) . Employing the standard assumptions (which will be made through- 
out this paper) of a 1.4 M Q , 10 km radius neutron star, LeiM = 2.7 x 10 38 ergs s^ 1 . In Figure^we plot the time evolution 
of L/LEdd over the first 1600 s of the 4U 1820-30 superbu rst. The observed luminosity was calculated by measuring the 0.1- 
40 keV flux predicted from the best fitting spectral models of Ballantvne & Strohmaver ( 2004), and assuming isotropic radiation. 
Separate curves are plotted for the diff erent distance e stimates described above. Since large super-Eddington luminosities are 
not expected for X-ray bursts (e.g., Joss & Melia 1987), this figure seems to indicate that the 5.8 kpc distance may be the most 
reasonable. However, by using equation^we find that increasing the mass and radius of the neutron star to, for example, 1.7 M Q 
and 12 km, would make the 6.4 kpc distance consistent with Lsdd- The 7.6 kpc distance is more challenging, but a much larger 
mass and radius (for example, M=2.3 Mq and i?=18 km) can allow the Eddington luminosity to approach the observed value. 
Such a large mass is implied by the highest frequency of the kHz QPO disco vered fro m 4U 1820-30 if it is assumed to corre- 
spond to the orbital frequency at the innermost stable circular orbit (Zhang et al. 1998). Finally, it is worth noting that since the 
superburst had a PRE phase to it, the emitting photosphere may be larger than the true stellar radius, resulting in LEdd being 
underestimated, and therefore overestimating L/L-^dd- Because of these uncertainties, we will, if necessary, use all three values 
of the distance in the following sections, and continue to assume iEdd = 2.7 x 10 38 ergs s _1 . 

X-ray bursts are seen from 4U 1820-30 when the persistent X-ray flux is < 4 x 10~ 9 ergs cm~ 2 s _1 , which is the low end 
of the observed range. This emission is powered by accretion onto the neutron star, so the flux can be converted to an accretion 
rate, M w 1.1 x 10 17 (d/6.4 kpc) 2 g s" 1 = 1.7 x 10~ 9 (d/6.4 kpc) 2 yr _1 (IStrohmaver & Browdl2002l) . consistent with 
the prediction of gravitational radiation driven mass-loss (Stel la et alJI 19871) . Using the above value of L E dd and the observed 
persistent luminosity, we obtain an Eddington fraction of M/Medd = 0.07(d/6.4 kpc) 2 , where M-Edd — -^Edd/^c 2 is the 
Eddington mass accretion rate (rj is the accretion efficiency). 



3. OUTFLOWS 
3.1. Radiatively driven winds 

Type I X-ray bursts and superbursts may sometimes reach the Eddington luminosity ^Edd, as witn essed in the radiation ex- 
pansion phase, where the photosphe re of the burning material can be lifted off of the neutron star (e.g.. lLewin. Vacca & Basinskal 
1984; Strohmaver & Bildsten 2003). At these luminosities the momentum in the radiation field may also affect the inner regions 
of the accretion disk, possibly even blowing it away from the central object. Below, we check if a radiatively driven wind can 
explain the behavior of 4U 1820-30 during the superburst. 

3.1.1. Estimates from Accretion Disk Theory 

Before d iving into the models of ra diative winds, it is worth calculating some simple estimates of outflow rates and energetics 
from basic Sha kura & Sunvae'vl ( ll973l) accretion disk theory to compare against the computations. 

The X-ray modeling showe d that the reflecting region moved from ~ 20 r g to ~ 100 r g in 500-1000 s 
(Balla ntvne & Strohm aver 2003). It is unknown what occurred to the material within 20 r g to prevent it from reflecting (this 
is speculated upon in §|6j, so we will only consider the evolution of the disk between 20 and 100 r g . The accretion rate onto the 
neutron star is low enoug h (for all three values of the distance) that almost all of the radiation pressure dominated region of the 
disk will be within 20 r g ( Shakura & Sunvaevl ll973l) . thus we will work only with the equations for a gas supported disk with 
electron scattering opacity. In this domain, the surface density of the disk can be written as 

S = 2.5 x 10 13 (^m p ) 4 / 5 (0.2(l + X))^^ (JjL^j " A/ 3 / 5 (J?j " J(R) 3 / 5 g cm" 2 , (2) 

where R is the radius along the disk, /i is the mean molecular weight of the accreting gas, m p is the proton mass, X is the mass 
fraction of hydrogen, a is the Shakura-Sunyaev viscosity parameter, and J(R) = (1 — (6r g /i?) 1//2 ) ("e.g.. lSvensson & Zdziarskil 

1994). For4U 1820-30, M « 1.1 X 10 17 (d/6.4 kpc) 2 g s~\ X « and/x w 1.3 for a fully ionized gas with acosmic abundance 
of heavy metals. This results in 

S . . x (J£r) ~"" ( - * ) " 5 (-?-) "" (*) ~"* J(Rf" g (3) 

\1AM @ J \ 1.1 x 10 17 gs-M V 6 - 4k P c / \ r gJ 

We then integrate eq.|3]between 20 r g and 100 r g to estimate the total mass in the disk that needs to be removed by an outflow, 
m = 3.2 x 10 19 (rf/6.4 kpc) 6 / 5 g, where we have assumed a = 0.1. If this mass is removed in ~ 500 s, then we obtain a mean 
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mass outflow rate of m = 6.4 x 10 16 (<i/6.4 kpc) 6 / 5 g s _1 . This is the same order of magnitude as the accretion rate onto the 
neutron star, which implies that an even higher outflow rate is needed in order to remove the inner disk. 

If this material is initially being lifted off an area A of the disk at the sound speed c s then by mass continuity, m = pc s A = 
nrrift e c s A where n is the mean number density of helium at the disk surface. Writing the soun d speed as c = \/kT I um n , and 
taking kT = 2.6 keV, the temperature of the blackbody near the beginning of the superburst (Strohmaver & Brown 2002), gives 
c s ~ 435 km s _1 . If the launching area is an ann ulus between 20 and 100 r„ , then n = 1.7 x 10 17 (d/6.4 kpc) 6 / 5 cm -3 . This is 
much lower than the number density predicted by Sha kura & Sunvaevl ill 9731) . but may be reasonable for the surface layers where 
the density must fall off rapidly. 

Equation|5]can also be used to calculate the gravitational potential energy of the disk material between 20 and 100 r g , 5.2 x 
10 38 ((i/6.4 kpc) 16 / 5 ergs. The energy released during the first ~ 1500 s of the 4U 1820-30 superburst is 10 41 ~ 42 ergs, assuming 
isotropic radiation. Therefore, there seems to be enough energy in the radiation field to completely remove the reflecting region 
of the accretion disk in a radiatively driven wind. The viability of such a wind is considered in more detail below. 

3.1.2. Numerical Models 

A numerical wind model can be used to calculate the radiative acceleration, velocity, and escape time for outflowing gas. For 
this paper, we define a model of a continuum-driven wind, launched from the accretion disk, and moving along radial streamlines 
(radial because the wind is driven by momentum taken from radially-streaming photons). This outflow structure therefore fills 
an open bi-cone in which gas flows radially outward, nearly along the surface of the accretion disk (for simplicity the wind is set 
here to flow at a constant 5° above the disk). 

The wind's acceleration is computed by taking the local velocity, density, ionization structure, and transmitted continuum into 
account. The outflow is launched from a user-defined radius with a nearly optically thick column at the disk (the maximum 
column that a radiation driving could support), with the initial outflow speed set to the sound speed. The photoionization of this 
outflowing gas is then simulated using Cloudy (Ferlan dEOQl at various heights in the wind. These simulations yield both the 
transmitted continuum and the ionization state of the gas from which the radiative acceleration can be calculated. Integrating this 
acceleration then yields the solution for the velocity as a function of distance along the wind. We then iterate on this solution to 
ensure that the density, photoionization, acceleration, and velocity profiles are all consistent 1 . 

The initial c onditions f or the model are chosen to match the observed and modeled properties of 4U 1820-30 
(Strohmaver & Bro wnl2002t iBallantvne & Strohmave d2004l). Fir st, the gas is illuminated with a 2.6 keV blackbody continuum, 
matching the temperature inferred by Strohmaver & Brown ( 2002) near the peak of the observed flux. We use gas abundances for 
cool extreme helium stars as given by Pande v et al.l (120011) . The initial number density was set to 10 17 cm~ 3 at the disk's surface, 
as implied by the simple calculation in the previous section 2 . The innermost streamline of the wind is set to a radius of 20 r g , 
the initial distance derived from the Fe Ka line. As mentioned above, the base of the wind is simulated as being approximately 
optically thick with Nu = 10 24 cm -2 to model the entire wind that would be launched radiatively from the disk and to set a 
conservative limit on the velocity of the outflowing gas. At the specified initial density, this column implies wind launching radii 
of 20 to ^70 r g , approximately the same range as the observed variation in the Fe Ka emission radius. These models are then 
run for a range of L/L^dd values. For these calculations, it is important to note that the LEdd defined in §|2]is replaced by the 
local Eddington luminosity for the accretion disk gas at r > 50r g where relativistic effects are not important. At those distances, 
we can effectively disregard the final bracketed expression in eq.^ yielding ZvEdd = 3-5 x 10 38 ergs s~ 4 . Thus, L/L^dd is set 
to the local value that the disk and wind would see. 

With the above parameters, models are then run for a range of Eddington ratios to examine the resultant outflow velocities. 
These simulations show that the gas is highly ionized by the superburst, such that the wind is accelerated only by electron 
scattering and bound-free opacity of the wind, justifying a posteriori the approximation that the wind is continuum-driven. 
Figure|2Jshows the velocity of this wind at ~ 100 r g as a function of the Eddington ratio. 

The time required for the gas to be launched from 20 to 100 r g is therefore calculated by a numerical integral over dr/v(s) 
where v(s) is the velocity as a function of distance along the outflow streamline. The results are shown in Figure[3]for the model 
where Nu = 10 24 cm~ 2 at the base of the wind. The gas would easily be driven beyond 100 r g in 500-1000 s as long as 
the source has h/L^dd ^ 1. These results are insensitive to changes in initial conditions, such as the initial density or blackbody 
temperature, since variations in either do not substantially change the ionization state of the gas or the continuum acceleration. 

We now ask whether the mass outflow rate from the disk is enough to remove a significant amount of mass from the accretion 
disk. From the above model with riHc = 10 17 cm -3 , for L/L-^dd ~ 1, u ~ 0.11c and n W j n d ~ 1.9 x 10 14 cm~ 3 at r = 100 r g = 
2.0 7 x 10 7 cm. Calculating m out = Airr 2 f pv where / is the covering fraction of the outflow, we find m out ~ 2.3 x 10 16 /gs _1 . 
In § 13.1.11 we estimated from standard accretion theory that an outflow rate of ~ 6 x 10 16 g s~ x would be required to remove the 
inner part of the disk. The geometry of the calculated streamlines imply / ~ 0.1, resulting in an outflow rate ~ 30 times smaller 
than the theoretical estimate. Therefore, these calculations show that it is feasible for a continuum-driven outflow to remove a 
significant portion of the inner accretion disk during a superburst, but only for L/L-Edd ^ 1 and it would not completely empty 
the inner accretion disk around 4U 1820-30 on the timescales observed (~ 500 s). The observed mass accretion rate for this 
system is 1.1 x 10 17 g s -1 (<i/6.4 kpc) 2 (§0, implying that the replenishment rate in the disk is very competitive with the outflow 

1 This model follows largely the same method as presented in Everett & Ballantvne 12004). In contrast to that earlier work, however, the angular momentum 
of the orbiting gas in these simulations is set solely by the mass of the central object (since the initial angular momentum was set by the quiescent luminosity, not 
the outburst luminosity modeled here). 

2 Our assumed base density is much greater than the 10 13 cm~ 3 limit beyond which Cloudy's three-body recombination calculations are defined fferland 
2002). However, for these luminosities, the ionization level is so high (U = ni on /njjo ~ 10 4 ' 5 , where n\ OYL is the density of ionizing photons) that this will not 
substantially affect our results. 
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rate. 



3.2. Thermally driven winds 

The gas near the surface of the accretion disk around the neutron star will be subject to intense X-ray heating by the explosion. 
If the temperature of the material increases to the point where the local sound speed exceeds the escape velocity then a thermal, 
pressure-driven wind could form, removing materia l from the disk. T his mechanism has been proposed as a means to disperse 
circumstellar gas around newly forming stars (e.g.. lFont et al]l2004l) . However, it is very unlikely that such a wind will form 
close to a compact object such as a neutron star, where the escape velocity will be a large fraction of the speed of light. This 
can be seen explicitly by comparing the temperature of the irradiated gas to the virial temperature kT v ~ GM fj,m p /R. During 
an outburst by the central source, the temperature at the accretion disk surface will be dominated by photoionization, so will 
maximize at the Compton temperatur e. For the case of 4U 1820-30, the peak temperature of the illuminating blackbody was 
~ 3 keV ( Strohmaver & Brown 2002), which will also be the value of the Compton temperature of the illuminated gas (§|j). For 
gas with /i w 1.3 orbiting a 1.4 Mq neutron star at 20 r g , kT v ~6x 10 4 keV. Clearly, thermal pressure alone will not cause the 
heated ga s to flow from the system. Other forces are required to overcome the strong gravity of the neutron star and cause gas to 

outflow ( $n\ . 

4. INFLOW 

X-ray photons produced in the explos ion on the surface of the neutron star can rem ove angular momentum from the accretion 
flow by Poynting-Robertson drag (Povnting 1903; Roberts on! 1 937t iBlumenthal 1974). This will increase the local accretion rate 
and potentially could drain mater ial from the inner part of the disk inward and onto the star JWalker & Meszarosl ll989t [Walker 
ll992UMiller & Lambll993l ll996). If this effect operates efficiently, it may provide another possible explanation for the apparent 
accretion disk evolution observed in 4U 1820-30. 

A simple estimate of the mass affected by radiation drag is when it interacts with its own rest mass of light, Mpp ~ Lt/c 2 , 

ilHD 



where t is the duration of the event (e.g., Walker & Meszaros 1989). In the case of 4U 1820-30, Lt ~ 10 4 



erg (§ 



which gives Mpr ~ (1.1 x 10 20 ~ 21 )/ g, where / is the fraction of the luminos ity inte rcepted by the disk. Even if / ~ 0.1, this 
estimate is comparable to the mass in the affected region of the accretion disk (S I3.1.H . We now check the timescale for angular 
momentum loss. 

For circular orbits with velocity v^, the rate of change of angular momentum, I, is related to the change of radius r by 



1 dl 
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(e.g.. lBlumenthalll97 



■Robertson drag gives the angular momentum loss rate as a function of the source luminosity L and r 
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where op is the Thomson cross-section, /3 = v$ /c, £ is the mean number of scattering el ectrons per atom ic mass units in the gas, 
/ w 5 is the increase in the angular momentum loss rate due to relativistic effects jMiller & Lambl 19961) . and 3 is a function of 
photon energy and contains corrections for Coulomb scattering (Blumenthal 1974). Substituting eq.[5]into eq.|4]and integrating 
gives the time for matter to spiral from a radius rg to the stellar surface at : 
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(6) 



where m is the mass of the scattering particle in the accretion flow. 

During t he 4U 1820-30 superburst the illuminating spectrum was a blackbody with a peak temperature of kT w 2.9 keV 
JStrohmaver & Brownl2002|). Since the vast m ajority of the illuminating photons have energies far above the ionization energy of 
helium (54.4 eV) ( — 0.5. Blumenthal (119741) shows in his Equation (29) a limiting expression for E in the case of a blackbody 
spectrum with kT <C mc 2 . Using this equation with kT = 2.9 keV gives S = 0.91. The radius expansion phase observed early 
on in the superburst implies the peak luminosity was ~ ^Edd, which we will take to be 2.7 x 10 38 ergs s _1 (§|2jl. Collecting all 
the above results and scaling the disk radii to gravitational radii r g gives a new expression for the time for matter to drain onto 
the neutron star due to Poynting-Robertson drag: 



t PR = 2 x 10" 
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(7) 



Figure|4]plots the predicted time for matter to spiral down from 100 r g to the stellar surface at ~ 5 r g during the first ~ 1600 s 
of the 4U 1820-30 superburst. The values of L/L^dd were taken from Fig.[0 and we show results for the different distance 
estimates to 4U 1820-30. For any distance, the Poynting-Robertson timescale is very small compared to the results of the 
spectral modeling which showed an evolution in the reflecting region over 500-1000 s. However, equation ignores processes 
such as outflows pushed by the radiation field (§ 13. 1> , viscous coupling of the gas, and radiative transfer effects in the optically 
thick disk material which will randomize the radiation field. All of these will slow down the infall rate onto the star, but it is 
unclear if the timescale can be increased by the four orders of magnitude needed to bring it in line with the observed behavior. 
Poynting-Robertson drag also suffers from the same dependence on the incident flux as the radiative winds in the previous section. 
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In contrast, the X-ray reflection results show th at the disk response was not dominated by the illuminating flux, but rather by the 
spectral shape (Ballantvne & Strohmaver 2004). 

Nevertheless, ther e are indications in the data that the Poynting-Robertson effect may be ongoing and affecting the disk struc- 
ture. iBallantvne & Strohmaverl (|2004) found that late in the burst decay phase, the inner radius of the disk (as judged by the 
reflection features) was relatively steady between 10 and 20 r g , much farther out than the inner-most stable circular orbit of 
~ 6 r g . At this point, the illuminating spectrum was cool and the observed flux was about 10 times lower than at the peak. It 
may be that at these low flux values, Poynting-Robertson is increasing the accretion rate at radii < 10 r g , thereby lowering the 
surface density and suppressing reflection features. 

5. LOCAL CHANGES IN DISK STRUCTURE 

The large X-ray flux released during a superburst or Type I X-ray burst will significantly heat the surface of the surrounding 
accretion disk. In this section we consider how the disk may respond to this heating 3 . 

For a Thomson-thin plasma dominated by photoionization, the maximum equilibrium temperature obtainable (although free- 
free cooling will somewhat reduce the maximum temperature) is the Compton temperature kTc, 

AkT c / u E dE = I u E EdE, (8) 

where Compton heating is balanced by Compton cooling. In the above expression is the spectral energy density of 
the radiation, and E is photon energy. For the case of a blackbody spectrum produced by a burst on a neutron star, 
Ue oc E 3 /(e E / kT - 1) and kT c « kT. In response to this increased temperature, the surface layers of the disk will have 
to adjust in order to maintain hydrostatic balance. The timescale of this readjustment is of the order of the dynamical time 
t z ~ t<j> ~ 4 x lO _6 (M/M0)(i?/r s ) 3//2 s (tFrank. Kin g & Raina l2002l) . which is essentially instantaneous at distances close to 
the neutron star. The new hydrostatic configuration will evolve to a larger disk scale height H oc c s r 3 / 2 oc T 1 / 2 r 3 / 2 , implying 
that the disk thickness will closely track the hardness of the burst spectrum. This is interesting because th e 4U 1820-30 X-ray 
reflec tion results show a correlation with the observed kT of the incident power-law. Using standard theory ( Shakura & Sunvaev 
119731) . the maximum accretion disk temperature for 4U 1820 -30 is 0.51 (d/6.4 kpc) -1 keV. The temperature of the blackbody 
spectrum peaked at ~ 2.9 keV ( Strohmav er & Brownll2002l) . so that at this point the scale-height of the disk atmosphere may 
have increased by a factor of at least 2. 

Significant changes to the surface density du e to internal viscos ity mechanisms occur on the viscous timescale, t V i SC , which is 
typically much lon ger tha n the dynamical time (Frank et al. 2002). For the gas pressure supported, electron scattering dominated 
disk described in ^ 13. 1.11 

t visc ~ 4.9 x 10 24 (0.2(1 + X))- 1 /5( iUmp )4/5 a -4/ 5 (^-) 8/5 M -2 / 5 (j?) J{R)-^ s, (9) 
where the symbols are the same as in eq.|3] Substituting values appropriate for4U 1820-30 (see §l2l& l3.14l l. we obtain 

Figure |5]plots i v ; sc as a function of radius, assuming three different values of the viscosity parameter a and a distance of 6.4 kpc 
(the other distances do not greatly affect the results). The viscous timescales predicted for a ~ 0.1-0.2 are in the appropriate 
range to potentially explain the evolution of the X-ray reflection features between 20 and 100 r„. The effective a can of course 
vary from reg ion to region in the disk and from time to time llHawlev & Krolik|l2001t lArmitage. Reynolds & Chiangll200ll 
Hawlev 2001), so this should be considered an average value. 

The viscous time can also be written as i v ; sc ~ a~ 1 (H/R)~ 2 t ( j > , which illustrates an inverse dependence on the disk scale 
height. The estimates of t v - lsc from equation[K)]assumed an unaltered scale height. As noted above, X-ray heating could increase 
H by ~ 2, thereby reducing the viscous time by ~ 4. If a is now assumed to be ~ 0.025-0.05 then this would maintain roughly 
the same i v ; sc as those plotted in Fig. [5] As t V i SC is changing as the gas is being heated, it may be difficult to obtain an accurate 
estimate of a. 

Since equation^lpredicts approximately the correct timescale for the observed evolution in 4U 1820-30, it strongly suggests 
that the superburst caused a significant change in the disk structure which is then propagated on a viscous time. X-ray heating 
will cause the disk scale height to increase, but this will occur locally very quickly. We hypothesize that the sudden, intense and 
relatively long-lasting X-ray illumination from a superburst can affect the angular momentum transport mechanism of the disk. 

6. DISCUSSION 

6.1. The Response of the Accretion Disk around 4U 1820-30 

This paper has discussed mechanisms that could influence the structure or dynamics of an accretion disk around a neutron star 
undergoing a Type I X-ray burst or a superburst. The moti vation was to determine an expla nation for the apparent evolution of 
the accretion disk around 4U 1820-30 during a superburst (Ballantvne & Strohmaver 2004). Here, we will look at the different 
possibilities identified above and determine their effectiveness in explaining the observations. 



3 This was discussed briefly by Ballantvne & Strohmaver 1 2004), but is presented more fully here. Values derived here supersede the earlier ones. 
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A radiative driven wind off the accretion disk is an effective means to remove mass from the inner regions of the disk. However, 
this scenario faces several difficulties when confronted with the data. First, the outflow can only be driven by super-Eddington 
luminosities. As mentioned in §|2| X-ray bursts are not expected to greatly exceed Lkhh, and so the superburst woul d have only 
reached that luminosity for a very short time at the initial stage of the explosion. Ballantvne & Strohmaver (2004) found that 
the X-ray reflecting region drastically changed ~ 500 s into the burst, when the luminosity would have been sub-Eddington. 
Second, the mass outflow rates predicted by the numerical wind calculations are too small to remove the inner disk (assuming 
it is a Shukura-Sunyaev disk) unless it subtends a significant fraction of the sky as seen from the X-ray source. Furthermore, 
the accretion rate through the disk is of the same order or even greater than the predicted outflow rates, and thus increases the 
difficulty in removing the inner disk. 

On the other hand, it may not be necessary to completely excavate the inner disk in order to remove the reflection features. 
All that is required is to replace the optically-thick, efficiently reflecting disk with something that is a poor reflector. If a wind 
is blown off the disk, then it may obscure the surface from the X-rays and thereby destroy the reflection features. The longer 
observed timescale may be related to this more inefficient evaporation of the disk. More detailed simulations are required to test 
this idea. 

In §|4]we found that Poynting-Robertson drag is an efficient mechanism to remove angular momentum from the gas. However, 
the predicted infall timescale is so rapid that the superburst would drain the entire accretion disk onto the star in only a few 
seconds. Clearly, the simple analytic estimate employed in §E]is not applicable to the case of optically thick accretion disks. One 
important issue with the analysis is the single particle assumption in equation^ If the illuminated gas on the disk is significantly 
coupled by the viscous stresses (e.g., magnetically) to other areas in the disk, then the mass and cross-section of this larger, 
more massive parcel of gas would replace the single particle in eq. The details of this coupling would obviously depend on 
the viscous physics within the accretion disk, but could produce a mechanism to reduce the angular momentum loss rate in this 
material. A full radiation hydrodynamics simulation which included radiative transfer would be needed to completely elucidate 
the effects of radiation drag on the disk material. It seems likely, there fore, that r adiation drag will only be important for the 
very innermost regions of the accretion flow (as concentrated on by Miller & Lamb 1996). The high luminosity of the burst will 
very quickly reduce the angular momentum from this material and increase the local accretion rate. We would expect the decay 
timescale to increase dramatically beyond this small r ange of radii. This is a potential explanation for why the inner radius of the 
reflecting region never falls below ~ 10 r g (see also § 16 2l below). 

The final mechanism considered in the paper involved neither outflow nor inflow (although it seems likely that both are occur- 
ring at some level), but rather a change in disk structure or surface density due to X-ray heating. This scenario is the easiest to 
explain the correlations seen in the X-ray data, as it predicts a strong dependence on the blackbody temperature (rather than the 
flux), and the viscous timescale is of the right order (with a ~ 0.03-0.2) to account for the trends seen in the reflections fits. 
However, it is the hardest to justify theoretically because the X-ray heating would change the local disk structure on a dynamical 
time, much faster than the viscous time. We suggest therefore that the intense and long lasting heating does have a large impact 
on the disk transport mechanisms which will change on the appropriate timescale. If the surface density drops globally over the 
inner accretion disk, due to the material being 'puffed up' by X-ray heating, then it could easily become an inefficient reflector. 
Unfortunately, to truly determine if this mechanism is the correct explanation, and to determine its importance relative to the 
inflow and outflow scenarios described above, a full numerical simulation which encompasses all the relevant physics will need 
to be performed. This is planned for future work. 

Of c ourse, there is the possibility that the spectral features identified by Strohma yer & Brownl ( 120021) are not due to reflection 
at all. Shaposhnikov & Titarchuk (2002) argue that since the disk will be heated very quickly to the blackbody temperature by 
the burst, all the heavy metals will be entirely ionized and unable to produc e reflection features. Thes e authors propose that 
the Fe Ka line and edge are formed in an outflow from the stellar surface ( Titarchuk & Shrader 2004). While it is possible 
that spectral features may form in an outflow (although it would be unusual to find a large c olumn of heavy metals close to the 
surface of a neutron star; Bildsten. Salneter & Wasserman 119921: iBildsten. Chang & Paerell 2003) it is not clear that the disk 
will be entirely ionized by the superburst. Numerical simulations of accretion disks su bject to the magnetorota tional instability 
consistently find that the gas is inhomogeneous and clumpy due to turbulence (e.g., [H awlev & Krolik 2001). These density 
inhomogeneities may enh ance the line emission in the X-r ay reflection spectrum (Ballantvne. Turner & Blaes 2004). 

The spectral fitting by Ballantvne & Strohmaver (2004) measured the ionization parameter £ = 47rFx/riHe of the reflecting 
region on the disk. Shakura-Sunyaev disk theory can be used to estimate £ and determine if the disk is expected to be completely 
overionized. Writing 7iHc = S/niHc-ff and Fx = L/4ttR 2 , we obtain 

£ w (11) 

The scale-height H for a gas-pressure dominated disk with electron scattering opacity is 

/ M \ 7/10 / R\ 21/20 

H = 1.4 X 10- 10 (/im p )- 2 / 5 (0.2(l + X)) 1 / 1 ^- 1 / 10 L^-J (y) M^JiR) 1 / 5 cm. (12) 
Using this expression, equation[3] X sa and (i k 1.3, we now obtain 

C « 1500C/AO (_JLJ) - 9/1 ° (*) - 7/2 ° ( *_\ ( _JL_) - 4/5 m -W erg s cm ( 13 ) 

Figure |6]plots the predicted £ at the radius measured by the spectral fitting of Bal lantvne & Strohmaverl ( T2004I) . assuming it is 
illuminated by the luminosity plotted in Fig.^ Following the results of §|5j we have also assumed a = 0.1. Overplotted on the 
figure are the measured values of the ionization parameter as inferred from the spectral fitting. 
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In contrast to the expectations of Shanos hnikov & Titarchukl (120021) . the disk does not have to be overionized by the burst. 
Rather, Figure[6]shows that the predicted £ is generally much smaller than the measured one, especially near the beginning of the 
superburst. A simple way to bring the theoretical prediction of £ in line with the measured one is to have a variable surface density 
which is ~ 40 times smaller over the first ~ 1000 s, but then becomes only ~ 8 times smaller a couple of hours later, as the burst 
is decaying. Of course, increasing the scale height also results in a larger £. The exact behavior of how the surface density would 
evolve as a function of radius and burst luminosity would have to be calculated with a numerical simulation. Finally, Figure [6] 
shows that the predicted £ follows the same trend as those inferred from the spectral fitting, supporting the ionized reflection 
interpretation of the observed spectral features. 

There are additional dynamical effects that arise in bursting systems that could affect the accretion flow. For example, while 
models of PRE bursts show that the luminosity from the photosp here will barely exce ed the Eddington l imit, they also in dicate that 
winds will be driven from the surface of the star ( Ebisuzaki, Hanawa & Sugi motdl 19831 iKatoll 19831 119851 lOuinn & Paczvhskl 
119851 Paczvrisk i& Anderson! 1 98(1 IPaczvnski & Proszvriski 1986; Joss & Melia 1987). Thus, while the burst radiation may not 
blow away the disk, there may be a mechanical interaction between the stellar wind and the accretion flow. However, the outflows 
generally carry very little kinetic energy compared to the radiative energy released in the burst, so any mechanical interaction is 
unlikely to be significant. Secondly, the X-ray heated gas on the disk surface will want to move to larger radii due to the gas 
pressure gradient as well as radiation pressure. This tendency to rotate at larger radii will be counteracted by viscous transport 
within the hot gas and radiative drag. Therefore, there may be some interesting dynamics within the disk as a response to X-ray 
bursts. 

6.2. Disk Truncation and Magnetic Effects 

Modeling of the reflectio n features in the spectra of the 4 U 1820-30 superburst showed that there was little to no reflection 
within a radius of ~ 15 r g (Ballantvne & Strohmaver 2004), much larger than the innermost stable circular orbit at 6 r g . One 
explanation for this apparently truncated accretion flow is locally enhanced accretion due to Poynting-Robertson drag (§ |4jl. 
Alternatively, this radius may correspond to the magnetospheric radius, R m , where the magnetic pressure from the neutron star's 
magnetic field is equal to the gas ram pressure in the accretion flow. Within R m the dynamics of the accreting gas is dominated 
by the magnetic field, and therefore the disk may be easily disrupted into a configuration which is unlikely to produce reflection 
features. 

The magnetospheric radius is dGhosh & Lambll978t fl979) 

where fi m is the magnetic moment of the stellar field (assuming dipole geometry). Substituting M = 1.4 M Q and M w 1.1 x 
10 17 (d/6.4kpc) 2 gs _1 , and taking r m « 15 r g , we find that fj, m » 2.8 x 10 26 (rf/6.4 kpc) G cm 3 , or B - 3x10 s (d/6.4kpc)G. 
This is a reasonable field strength for a non-pulsating LMXB system dPsaltis & Ch akrabartv 1999). In fact, almost all systems 
which exhibit X-ray bursts, and LMXBs in general, do not show X-ray pulsations (Lewin et al. 1993), implying much weaker 
field strengths then pulsars (which have B ~ 10 12 G). Thus, it is possible that the lack of X-ray reflection from within ~ 15 r g 
may be due to magnetic truncation. 

6.3. Prospects for Future Observations 

The superburst from 4U 1820-30 occurred when the system was being observed by RXTE. As a result, the sensitive narrow- 
band instruments onboard were able to pick up the refl ection features and track their evolution as the burst ev olved. So far only 
one other system has been as fortunate (4U 1636-53: IStrohmaver & Markwardt 2002; Kuulker s' et alJl2004l) . Superbursts are 
the ideal events for which to search for X-ray reflection features from the accretion disk as they last for a number of hours. 
Moreover, because of this long timescale, the analyses presented in this paper predict that we are more likely to observe changes 
to the accretion disk during superbursts than during the much shorter Type I X-ray bursts. The major difficulty is the unknown 
recurrence time o f superbursts which can vary over a number of years depending on the composition of the accreted material 
dCumm ing 2003). Nevertheless, we suggest that rapid X-ray followup of superbursts would be an important secondary science 
objective for the Swift mission. Detecting reflection features from a variety of systems with different orbital and spin periods, 
magnetic field strengths and evolutionary paths would allow examination of the accretion process over a wide parameter space. 

As the overall geometry is the same and the light travel time is much shorter than the burst duration, X-ray reflection features 
should be produced in Type I bursts as well. The difficulty will be detecting them with sensitive enough instruments to do time 
averaged spectral analysis let alone time-resolved studies. However, Type I events are much more common than superbursts 
allowing repeatable experiments and tests of hypotheses. Sensitive spectral analysis of Type I X-ray bursts would be an important 
objective for a future high-throughput observatory. 

7. CONCLUSIONS 

Type I X-ray bursts or superbursts on neutron stars produce a significant amount of energy (~ 10 39 or ~ 10 42 ergs, respectively) 
in a short amount of time. In this paper, we have discussed how this energy, in the form of a X-ray flux, would affect the 
surroundin g accretion flow. The goa l was to elucidate the meaning behind the apparent evolution of the accretion disk uncovered 
by Ballantvne & Strohmaver (2004) in their analysis of X-ray reflection features observed during the 4U 1820-30 superburst. 

Three separate physical processes were considered. The first, radiatively driven winds, was found to be occur, but only for 
L/L-&dd ^ 1 and, for a thin accretion disc, the predicted outflow rate was less than the accretion rate, so this effect may not 
produce a large change to the accretion flow. Poynting-Robertson drag, the second process, was in fact predicted to be too 
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efficient, and would have removed the inner accretion disk is less than a second. This was the result of using a too simple 
description of the problem, but Poynting-Robertson may be dominant very close to the central star. The final process considered 
was a change in surface density due to X-ray heating. This successfully reproduced the observed timescale for a Shakura-Sunyaev 
viscosity parameter of a ~ 0.03-0.2, but requires the ansatz that the X-ray heating will result in significant changes to the surface 
density of the disk that would propagate on a viscous time. Our conclusion is that, while the X-ray heating may dominate the 
explanation of the observations, all three processes will be ongoing during a superburst, and numerical simulations are required 
to fully understand the response of the accretion disk to such large explosions. 

Finally, we emphasize that future rapid follow-up observations of X-ray bursts (both super and otherwise) will be vital in the 
progress of this problem. X-ray reflection features should be common in these systems and any observed changes are likely to 
correspond to changes in the accretion flow. Further examples are needed to test the ideas presented here as well as provide 
constraints on accretion theory. 

The authors thank C. Matzner and T. Strohmayer for useful discussions and comments, and the anonymous referee for sug- 
gestions which improved the presentation of the paper. This research was supported by the Natural Sciences and Engineering 
Research Council of Canada. 
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FIG. 1. — Time evolution of L/Le^cI during the first ~1700 s of the the 4U 1820-30 superburst. The three curves show the results for different distance 
estimates to the system. I/Edd is taken to be 2.7 X 10 38 ergs s — 1 , assuming a canonical neutron star (M = 1.4 Mp), R = 10 km) w ith a helium atmosphere. 
The observed luminosity was calculated between 0. l^M) keV using the best-fitting spectral models of Ballantvne & Strohmaver 1 2004 ) and assuming isotropic 
radiation. 
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FIG. 4. — Timescale for matter to drain onto the neutron star via Poynting-Robertson drag during the initial ~ 1600 s of the 4U 1820-30 superburst. As 
the timescale is directly related to the luminosity, three lines are plotted for different estimates of the distance to the system. The luminosity was calculated as 
described in the caption to Fig.Q In all three cases, t-p-R, is much less than 1 s. 



Suddenly Heated Accretion Disks around Neutron Stars 



13 




20 40 60 80 100 

Radius (r ) 

FIG. 5. — Viscous timescale for a gas pressure supported disk with electron scattering opacity as a function of disk radius. The three curves denote different 
values of the Shakura-Sunyaev viscosity parameter. A d istan ce to 4U 1 820-30 of 6.4 kpc was assumed, although use of the other two distances do not significantly 
change the results. The curves were calculated with Eg. 1101 and thus assume an unaltered scale height H. 
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FIG. 6. — The dashed line plots the predicte d value of the disk ionization par ameter from equation ll3l during the 4U 1820-30 superburst. The points are the 
best fit values from the reflection modeling of Ballantvne & Strohmaver 1 2004). A distance of 6.4 kpc and a = 0.1 was assumed for the theoretical estimate. 
In order to bring the prediction in line with the observations, the surface density must be decreased by about a factor of 40 at the beginning of the burst, but this 
decreases to ~ 8 at later times. 



